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ORIGINAL ARTICLE 

Evaluation of 4p-Hydroxycholesterol as a Clinical 
Biomarker of CYP3A4 Drug Interactions Using a Bayesian 
Mechanism-Based Pharmacometric Model 

TA Leil 1 , S Kasichayanula 1 , DW Boulton 1 and F LaCreta 1 

A Bayesian mechanism-based pharmacokinetic/pharmacodynamic model of cytochrome P450 3A4 (CYP3A4) activity 
was developed based on a clinical study of the effects of ketoconazole and rifampin on midazolam exposure and plasma 
4|3-hydroxycholesterol (4|$HC) concentrations. Simulations from the model demonstrated that the dynamic range of 4|}HC as a 
biomarker of CYP3A4 induction or inhibition was narrower than that of midazolam; an inhibitor that increases midazolam area 
under the curve by 20-fold may only result in a 20% decrease in 4|}HC after 14 days of dosing. Likewise, an inducer that elevates 
CYP3A4 activity by 1.2-fold would reduce the area under the curve of midazolam by 50% but would only increase 4|3HC levels 
by 20% after 14 days of dosing. Elevation in 40HC could be reliably detected with a twofold induction in CYP3A4 activity with 
study sample sizes (N ~ 6-20) typically used in early clinical development. Only a strong CYP3A4 inhibitor (e.g., ketoconazole) 
could be detected with similar sample sizes. 
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It is important to understand the potential of a new chemi- 
cal entity to have significant pharmacokinetic (PK) interac- 
tions with commonly administered comedications as early in 
the exploratory drug development process as possible. This 
knowledge enables important decisions to be made about 
either compound selection or revisions to the clinical develop- 
ment program before the investment of significant resources 
in intermediate- and late-phase clinical trials. PK interactions 
involving metabolic routes of drug elimination are generally 
well understood and tend to be the most straightforward to 
screen for using in vitro and animal models before selection 
of a compound for exploratory clinical development. Strong 
inhibitors and inducers of cytochrome P450 3A4 (CYP3A4) 
are generally excluded by in vitro and in vivo preclinical 
experiments before selection of a compound to move to first- 
in-human clinical studies. However, some compounds with 
the potential for clinically meaningful drug interactions do 
advance into human clinical studies, and some are eventu- 
ally approved for use in patients. 1 

The traditional approach for detecting the potential for 
CYP3A4-mediated drug interactions for compounds in 
clinical development is to conduct a clinical PK drug inter- 
action study in human subjects, using a probe substrate of 
CYP3A4, such as midazolam (MDZ). These studies are typi- 
cally required as part of the regulatory approval process for 
drugs in which clinically relevant CYP3A4 induction or inhibi- 
tion cannot be ruled out based on nonclinical experiments. 2 
4p-Hydroxycholesterol (4pHC) is an endogenous oxidized 
metabolite of cholesterol that is specifically produced via 
metabolism of cholesterol by CYP3A4 in the liver. 3 Induc- 
tion of CYP3A4 causes elevation in plasma concentrations 
of 4pHC, whereas inhibition of CYP3A4 leads to reduction. 45 
Therefore, 4pHC has the potential to be used in early-phase 
clinical trials as an endogenous biomarker to determine 



the potential for a compound to induce or inhibit CYP3A4. 
4oc-Hydroxycholesterol (4ocHC) is a cholesterol-oxidation 
product that is not produced through metabolism by CYP3A4. 
It has been proposed as a possible endogenous metabolite 
that could be used to normalize for differences in general 
metabolic activity between subjects, thus reducing the impact 
of intersubject variability on estimation of changes in 4pHC. 

Before the routine implementation of 4pHC as a biomarker 
in clinical development, it is important to understand its 
dynamic properties in response to well-accepted inducers 
and inhibitors of CYP3A4. In addition, the inter- and intra- 
subject variability in the steady-state levels of 4pHC and the 
response to CYP3A4 modulation in healthy subjects and 
patients will determine its reliability as a biomarker in clini- 
cal studies. Several experimental clinical studies in healthy 
subjects have already been conducted using strong inhibi- 
tors (atazanavir/ritonavir and itraconazole) 56 or inducers 
(rifampin, efavirenz, and carbamazepine) 7-9 of CYP3A4 to 
evaluate the potential of 4pHC to be used as a clinical bio- 
marker. The results of these studies suggest that because 
of the relatively long half-life of 4pHC (1-3 weeks) 10 and the 
fact that inducers result in elevation of 4pHC concentrations, 
it may be a more valuable measure of CYP3A4 induction as 
opposed to inhibition. Indeed, Yang and Rodrigues 11 demon- 
strated that these results would be expected based on a the- 
oretical model of 4pHC dynamics. However, a model-based 
characterization of the pharmacodynamics (PD) of 4pHC, 
and its variability, in humans will provide a thorough under- 
standing of its dynamic properties as a biomarker. A model 
developed in such a way would permit prospective prediction 
of the 4pHC biomarker response in clinical studies based on 
preclinical data. This would enable an early decision to be 
made regarding whether to collect this biomarker, as well as 
the timing and frequency of collection. 
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In this article, we report the development and utility of a 
Bayesian mechanism-based PK/PD model of 4pHC using 
data from a clinical study that evaluated 4pHC as a biomarker 
of both CYP3A4 induction by rifampin and CYP3A4 inhibition 
by ketoconazole. The objectives for developing this model 
based on mechanistic information and clinical data are as fol- 
lows: (i) to provide a tool that permits integration of nonclini- 
cal information on the potential of an new molecular entity 
to have CYP3A4 interaction, so that decisions can be made 
about how to use 4pHC in early clinical development; (ii) to 
provide a tool to conduct clinical trial simulations to optimize 
the design of clinical studies in which 4pHC is used to evalu- 
ate CYP3A4 inhibition/induction; and (iii) to evaluate the util- 
ity and dynamic range of 4pHC relative to MDZ as a measure 
of CYP3A4 induction/inhibition by compounds with different 
induction/inhibition potential. 



RESULTS 

Development of the PK/PD model for 4pHC 

A schematic diagram of the nonlinear mixed effects PK/ 
PD model for CYP3A4-mediated drug interactions is 
shown in Figure 1. The model describes the PK of MDZ 
and ketoconazole, the formation of 4pHC and 4ocHC from 
cholesterol (the PD measure), and the induction caused by 
600 mg of rifampin. 4ocHC is not a substrate for CYP3A4 
but may be useful in normalizing for differences in gen- 
eral metabolic rates among individuals when estimating 
changes in 4pHC. Separate compartments for induction 
and inhibition of CYP3A4 in the gut vs. the liver were 
implemented. Bayesian priors were used for all the model 
parameters. Informative priors were used, when possible, 
based on the availability of published data. Plasma keto- 
conazole concentration was not measured in the clinical 




Figure 1 Diagram of CYP3A4 semi-mechanistic PK/PD model for CYP3A4 interactions. 4aHC, 4a-hydroxycholesterol; 4pHC, 
4p-hydroxycholesterol; C Keto , concentration of ketoconazole; CL, clearance from central compartment; CL |NT , intrinsic clearance; C MDZ , 
concentration of MDZ; CYP, fold induction in CYP3A4; FPE, first-pass effect for MDZ; Fu, fraction unbound for MDZ; Fuk, fraction unbound for 
ketoconazole; /A/H 4J3HC , inhibition function for effect of ketoconazole on formation of 4p-hydroxycholesterol; INH MDZ , inhibition function for effect of 
ketoconazole on midazolam's intrinsic clearance; K, central compartment for ketoconazole; /C deg4 HC , elimination rate for 4a-hydroxycholesterol; 
K deg4(3HC , elimination rate for 4p-hydroxycholesterol; K, inhibition constant for ketoconazole; K m , concentration for half maximal metabolic 
clearance rate; K , peripheral compartment for ketoconazole; K , synthesis rate for 4p-hydroxycholesterol in the absence of induction or 
inhibition; K 3 4ocHC , synthesis rate for 4oc-hydroxycholesterol; K syn * synthesis rate for 4p-hydroxycholesterol; K c , synthesis rate for total 
cholesterol; K TR " gut transit rate constant; M c , central compartment for MDZ; M , peripheral compartment for MDZ; MTT, mean transit time in gut; 
N, number of compartments; Q H , hepatic blood flow; Rif, rifampin; V c , volume of central compartment; V mgx , maximum rate of metabolic clearance. 
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study; therefore, highly informative priors for ketoconazole 
PK parameters were used based on a previously published 
meta-analysis of ketoconazole PK (Table 1). Similarly, the 
priors for the PK parameters for MDZ and the interaction 
with ketoconazole were also based on a previously pub- 
lished meta-analysis (Table 1). The objective function and 
the fixed effects parameters appeared to converge to their 
optimal stationary values based on examination of the cat- 
erpillar plots and evaluation of the Gelman and Rubin inter- 
chain diagnostics for the three independent Markov chain 
Monte Carlo (MCMC) chains (Supplementary Figure S1 
and Supplementary Table S1). The prior and posterior 
means and SDs of the model parameters can be found 
in Table 1. Separate K. values were required to describe 
inhibition of CYP3A4 for MDZ (K = 0.125 nmol/l) vs. 4pHC 
(K = 50 nmol/l). Separate values were also required 
to describe induction of CYP3A4 by rifampin in the gut 
(29.6-fold) vs. the liver (3.43-fold). 



Evaluation of the PK/PD model for 4pHC 

The posterior predictions of MDZ plasma concentrations and 
4pHC serum concentrations from the model were compared 
with the observed data, both visually and by comparison of 
the summary statistical measures. The posterior distribution 
of predicted MDZ plasma concentrations vs. time compared 
with the observed data following oral and i.v. MDZ, either 
alone or in the presence of ketoconazole or rifampin, are 
shown in Figure 2a-f. The model is able to capture the time 
course of MDZ concentrations and the effects of CYP3A4 
inhibition by ketoconazole and induction by rifampin. Com- 
parison of the summary statistics for the posterior model pre- 
dictions of MDZ exposures to the observed statistics reveals 
some minor discrepancies in the model predictions (Table 2). 
There was a slight overprediction of the area under the curve 
(AUC) for oral MDZ in the presence of ketoconazole (-20%) 
and a slight overprediction of the AUC of oral MDZ in the 
presence of rifampin (-29%). However, the observed oral 



Table 1 Prior and posterior model parameters 



Parameter 


Description 


Prior mean (SD) 


Prior ISV (df) 


Pnctorinr 

mean (SD) 


Pnctofinr 

ISV (SE%) 


F M 


MDZ bioavailability 21 


0.69 (0.71) 


0.32(1) 


1.0 (0.1) 


0.42 (39) 


K Am 


MDZ absorption rate (1/h) 20 


1.86 (0.67) 


0.45(1) 


3.91 (0.42) 


0.61 (49) 


Q m 


MDZ intercompartmental clearance (l/h) 20 


42.1 (18) 


0.45 (4) 


13.4(1.43) 


0.52 (30) 


v Cm 


MDZ central volume (I) 20 


78.8 (30) 


0.33 (4) 


45.3 (6.4) 


0.87 (28) 


v Pm 


MDZ peripheral volume (I) 20 


53.3 (23) 


0.45 (4) 


114(18.4) 


1 .6 (32) 


v max m 


MDZ V max (umol/h) 20 


5,125.4(1,844) 


0.39 (40) 


1,100 (122) 


0.48 (21) 


K mm 


MDZ K m (umol/l) 20 


2.11 (0.06) 


0.032 (40) 


2.10 (0.064) 


0.037 (27) 


F us 


MDZ fraction unbound 21 


0.04 (0.001) 


0.032 (40) 


0.040 (0.0009) 


0.037 (28) 


Q 


Hepatic blood flow (l/h) 21 


80 (0.08) 


0.032 (40) 


80 (0.078) 


0.039 (30) 


F K 


Ketoconazole bioavailability 21 


0.70 (0.045) 


0.032 (40) 


0.70 (0.0002) 


0.037 (28) 


*a 


Ketoconazole absorption rate (1/h) 21 


0.61 (0.005) 


0.032 (40) 


0.61 (0.005) 


0.037 (28) 


Q 


Ketoconazole intercompartmental clearance (l/h) 21 


3.22 (0.1) 


0.032 (40) 


3.22 (0.1) 


0.037 (28) 




Ketoconazole central volume (I) 21 


20.6 (0.65) 


0.032 (40) 


20.4 (0.7) 


0.037 (28) 


K 


Ketoconazole peripheral volume (I) 21 


24.4 (0.77) 


0.032 (40) 


24.5 (0.79) 


0.037 (28) 


v max 


Ketoconazole V mgx (umol/h) 21 


23.1 (0.73) 


0.032 (40) 


23.0 (0.70) 


0.037 (28) 


K m 


Ketoconazole K m (umol/l) 20 


0.48 (0.004) 


0.032 (40) 


0.48 (0.004) 


0.037 (28) 




Ketoconazole F plasma 26 


0.029 (0.001) 


0.34 (4) 


0.029 (0.0007) 


0.032 (23) 


B_BHCH 


Baseline 4(3-hydroxycholesterol (nmol/l) 27 


29 (14) 


0.032 (40) 


81.2 (4.7) 


0.32 (28) 


BHCJHL 


Half-life of 4(3-hydroxycholesterol (h) 19 


234(118) 


0.032 (40) 


124 (7.6) 


0.20 (47) 


CH_HL 


Half-life of total cholesterol (h) 28 


192 (6.1) 


0.032 (40) 


192 (0) 


0.031 (1E-5) 


BCH 


Baseline cholesterol (umol/l) 27 


4,500 (2,240) 


0.031 (4) 


4,244(127) 


0.15(30) 


B_AHCH 


Baseline 4a-hydroxycholesterol (nmol/l) 27 


6.6 (3.3) 


0.42 (4) 


10.6 (0.53) 


0.22 (29) 


AHCJHL 


Half-life of 4a-hydroxycholesterol (h) 


234 (7.4) 


0.032 (40) 


234 (0) 


0.031 (1E-5) 


CYP_HL 


CYP3A4 half-life (h) 11 


70 (2.2) 


0.032 (40) 


70 (0) 


0.032 (1 E-5) 


Ki 


K. of CYP3A44 inhibition for MDZ (nmol/l) 21 


14(3.0) 


0.45(1) 


0.125 (0.226) 


2.2 (168) 


FOLDINDS 


Fold induction of systemic CYP3A44 11 


12(6) 


0.45(1) 


3.43 (0.46) 


0.42 (59) 


MTT 


Mean gut transit time for MDZ 


2.5 (1.6) 


0.45(1) 


0.544 (0.057) 


0.49 (34) 


FOLDINDG 


Fold induction of gut CYP3A4 


12(6) 


0.032 (40) 


29.6 (1.5) 


0.032 (24) 


K.4B 


K. of CYP3A4 inhibition for 4(3-hydroxycholesterol (nmol/l) 


14(3.0) 


0.45(1) 


50 (17) 


0.70 (130) 


G 1 


Proportional SD for MDZ 






0.36 (0.077) 




a 2 


Additive SD for MDZ (nmol/l) 






0.046 (0.037) 






Additive SD for cholesterol (nmol/l) 






1,170 (284) 




a 4 


Additive SD for 4(3-hydroxycholesterol (nmol/l) 






28.3 (7.2) 




°5 


Additive SD for 4a-hydroxycholesterol (nmol/l) 






4.19 (1.02) 





CYP3A4, cytochrome P450 3A4; df, degrees of freedom; ISV, intersubject variance; MDZ, midazolam. 
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Figure 2 Evaluation of posterior and prospective predictive performance of the PK/PD model. Observed and posterior predicted concentrations 
for MDZ and 4pHC for: (a) 2 mg oral MDZ; (b) 2 mg oral MDZ with ketoconazole; (c) 2 mg oral MDZ with rifampin; (d) 0.8 mg i.v. MDZ; (e) 0.8 mg 
i.v. MDZ with ketoconazole; (f) 0.8 mg i.v. MDZ with rifampin; (g) 4|3HC after 14 days of 400 mg q.d. ketoconazole; (h) 4PHC after 14 days of 
600 mg q.d. rifampin; (i) observed and prospectively predicted MDZ ratio for 13 compounds; and (j) observed and prospectively predicted 
4(3HC vs. time following carbamazepine administration. Plots a through h show the observed concentrations vs. time, the median and 90% 
intervals for the posterior predictions, and the 90% credible intervals. Plot i shows the observed vs. predicted with the line of unity (solid line) 
and the lines representing twofold prediction error (dotted lines). Plot j shows the pooled 4pHC observations from Wide etal., as well as the 
median and 95% interval of the population predictions. 4pHC, 4p-hydroxycholesterol; CYP3A4, cytochrome P450 3A4; MDZ, midazolam; PD, 
pharmacodynamics; PK, pharmacokinetics. 



and i.v. MDZ AUC ratios following ketoconazole or rifampin 
administration are within the credible intervals of the poste- 
rior predictions. 



The posterior predictions of 4pHC concentration vs. time 
compared with the observed data are shown in Figure 2g,h. 
The predictions of 4pHC levels following inhibition by 
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Table 2 Observed and posterior predicted midazolam exposures and 4(3-hydroxycholesterol 



Parameter 


Observed 3 (SD) 


Predicted 3 (95% CI) 


Oral midazolam AUC 0 _ 24 (nmol/l*h) 


72.8 (5.21) 


76.5 (62.3-96.1) 


Intravenous midazolam AUC 0 _ 24 (nmol/l*h) 


55.95 (10) 


61.1 (55.2-67.5) 


Oral midazolam AUC 0 24 (nmol/l*h) with ketoconazole 


963.7 (52.6) 


932 (805-1,080) 


Intravenous midazolam AUC 0 24 (nmol/l*h) with ketoconazole 


271 (31.9) 


339 (299-387) 


Oral midazolam AUC 0 _ 24 (nmol/l*h) with rifampin 


4.35 (0.619) 


6.09 (4.72-7.71) 


Intravenous midazolam AUC 0 _ 24 (nmol/l*h) with rifampin 


24.49 (2.29) 


25 (22.4-27.8) 


Oral midazolam ratio: ketoconazole 


11.41 (1.4) 


12.3 (9.13-15.9) 


Intravenous midazolam ratio: ketoconazole 


4.612 (0.538) 


4.84 (4.1-5.67) 


Oral midazolam ratio: rifampin 


0.06096 (0.00873) 


0.102 (0.0858-0.122) 


Intravenous midazolam ratio: rifampin 


0.4255 (0.0879) 


0.41 1 (0.37-0.448) 


Baseline 4(3-hydroxycholesterol 


81.9 (5.73) 


79.1 (69.9-87.5) 


Change from baseline in 4(3- hydroxy cholesterol with rifampin at day 14 


291 (24.5) 


214(166-314) 


4(3-Hydroxycholesterol ratio with rifampin 


3.91 (0.245) 


3.14(2.66-4.14) 


Change from baseline in 4(3-hydroxycholesterol with ketoconazole at day 14 


-32.5 (5.77) 


-27.7 (-31.9 to -23.5) 


4(3-Hydroxycholesterol ratio with ketoconazole 


0.675 (0.0577) 


0.723 (0.681-0.765) 



AUC, area under the curve; CI, confidence interval. 

a Summary statistics are geometric means for AUCs and 4(3-hydroxycholesterol change from baseline, and medians for the ratios. 



ketoconazole or induction by rifampin appear to be adequate. 
Table 2 provides the summary statistics of predicted and 
observed means for 4pHC serum concentrations at baseline 
and on day 14. 

To evaluate the prospective predictive performance of 
the model, predictions were generated for 13 known inhibi- 
tors and inducers of CYP3A4 that were not included in the 
analysis dataset. Figure 2i shows the observed vs. predicted 
MDZ AUC ratio for these compounds. The data and refer- 
ences used in the prospective predictions can be found in 
Supplementary Table S2. Figure 2i shows the predicted 
and observed time course for the change from baseline in 
4pHC following administration of carbamazepine (CYP3A4 
inducer) for 24 weeks. The observed 4pHC vs. time data for 
carbamazepine treatment come from a clinical study by Wide 
et a/. 9 in pediatric epilepsy patients. 

Evaluation of dynamic range of 4pHC vs. MDZ 

Simulations using the model were conducted to characterize 
the dynamic range of 4pHC as a biomarker of CYP3A4 induc- 
tion and inhibition in comparison to MDZ. Figure 3a shows the 
percent increase in 4pHC compared with the fold decrease in 
MDZ AUC in the presence of inducers with different CYP3A4 
induction potentials after 14 days of treatment. 4pHC does not 
appear to be as sensitive as MDZ for inducers; for an inducer 
with 1.05-fold CYP3A4 induction potential (equivalent to a 
105% increase in enzyme activity), the AUC of MDZ would be 
expected to be reduced by 20%, whereas the AUC of 4pHC 
would be expected to increase by only 4%. To see an increase 
in 4pHC of 20%, CYP3A4 induction would need to be 1.2- 
fold (equivalent to a 120% increase in enzyme activity), suf- 
ficient to reduce MDZ AUC by 50%. Figure 3b,c show the fold 
increase in MDZ AUC vs. the percent decrease in 4pHC in the 
presence of CYP3A4 inhibitors, with different K. values and a 
half-life of 2 h, administered as 400 mg q.d. or 200 mg b.i.d. The 
difference in dynamic range between 4pHC and MDZ is more 
pronounced for inhibitors than that for inducers. Inhibitors that 
reduce 4pHC by 20% would be expected to increase MDZ 
AUC by 20- to 26-fold, depending on the dosing regimen. 



Clinical trial simulations of CYP3A4 induction and 
inhibition 

To characterize the potential utility of 4pHC as a biomarker 
of CYP3A4 induction in clinical trials, simulations were con- 
ducted using compounds with different CYP3A4 induction 
potentials. Figure 4a shows the change in 4pHC vs. time 
for four CYP3A4 inducers with different CYP3A4 induction 
potentials, following 14 days of treatment. Induction potential 
is defined as the fold increase in CYP3A4 by day 1 4. Induction 
potentials of 1 .5-, 2-, 3-, and 4.6-fold would be expected to 
result in changes in 4pHC of 40, 75, 150, and 244%, respec- 
tively, on day 14. Simulation of clinical trials with different 
sample sizes show the relationship between sample size and 
the probability to detect a significant difference in 4pHC, rela- 
tive to placebo-treated subjects, for compounds with different 
induction potentials (Figure 4b). A significant difference was 
defined as lack of overlap between the 95% confidence inter- 
vals of the treated group relative to placebo, with the assump- 
tion of equal sample sizes between the groups. There is a 
relatively linear relationship between induction potential and 
the sample size to achieve 80% probability (Figure 4c). 

The change in 4pHC as a result of liver CYP3A4 inhibition 
in the model is driven by unbound concentrations of inhibitor 
in the systemic circulation as a function of time and there- 
fore depends on the K. for the enzyme, the PK properties, 
and the dosing regimen of the compound. Figure 5a shows 
the change in 4pHC vs. time for five CYP3A4 inhibitors with 
different PK properties (defined primarily by half-life), K. val- 
ues, and dosing regimens, following 14 days of treatment. 
Simulation of clinical trials with different sample sizes show 
the relationship between sample size and probability to 
detect a significant difference in 4pHC relative to placebo- 
treated subjects for CYP3A4 inhibitors with different proper- 
ties (Figure 5b). Because several compound-specific factors 
influence the change in 4pHC on day 14, an examination of 
different compound properties was conducted to determine 
the measure that would best correlate with the effect on 
4pHC. The concentration of inhibitor divided by the K. (I/K) 
appears to be the most useful parameter for predicting the 
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Figure 3 Comparison of dynamic range of 4pHC vs. MDZ after 14 days of treatment with inducers or inhibitors of CYP3A4. (a) Fold change 
in 4pHC vs. MDZ AUC ratio for different levels of CYP3A4 induction; blue dotted lines show an inducer with 1 .05-fold induction potential (0.8 
MDZ AUC ratio and 4% increase in 4pHC). Green dotted lines show an inducer with 1 .2-fold induction potential (0.5 MDZ AUC ratio and 20% 
increase in 4pHC); black dotted lines show an inducer with 4.6-fold induction potential (0.1 MDZ AUC ratio and 260% increase in 4PHC). 
(b) Fold change in 4pHC vs. MDZ AUC ratio for inhibitors with different CYP3A4 K. values and a 2-h half-life administered at 400 mg q.d. Red 
dotted lines show an inhibitor with an [f\/K. iCmin of 0.003 (20-fold increase in MDZ AUC ratio and 20% decrease in 4PHC). (c) Fold change in 
4pHC vs. MDZ AUC ratio for inhibitors with' different CYP3A4 K. values and a 2-h half-life administered at 200 mg b.i.d. Red dotted lines show 
an inhibitor with an [l\/K [Cm]n of 0.0012 (26-fold increase in MDZ AUC ratio and 20% decrease in 4PHC). 4|3HC, 4p-hydroxycholesterol; AUC, 
area under the curve; CYP3A4, cytochrome P450 3A4; MDZ, midazolam. 



change in 4pHC at the end of treatment. Both UK for aver- 
age plasma concentration at steady state (C avg ) and steady- 
state C min show similarly strong correlations with the change 
in 4pHC (data not shown). Figure 5c shows the relationship 
between UK for C min at steady state vs. the sample size to 
achieve 80% probability to detect a significant difference 
in 4pHC vs. placebo treatment. In contrast to the inducers, 
there is greater variability in this relationship because of the 
several factors involved in producing changes in 4pHC (i.e., 
CYP3A4 K, PK properties, and dosing regimen). It is impor- 
tant to determine that the type 1 error rate (i.e., false-positive 
rate) is consistent with expectations when interpreting the 
results of clinical trial simulations. Clinical trial simulations 
using the PK/PD model and assuming no induction/inhibition 
of CYP3A4 (n = 8 each in placebo and treated) determined 
the type 1 error to be -5%. 

DISCUSSION 

We have reported the development of a nonlinear mixed 
effects PK/PD model to evaluate the utility of 4pHC as a 



biomarker of induction and inhibition of CYP3A4 in clinical 
studies. A Bayesian approach was utilized for estimation of 
the model parameters, both to leverage prior knowledge in 
estimation of parameters based on the current clinical trial 
dataset and to compensate for data that was not collected in 
the study. The Bayesian approach is useful in facilitating the 
development of mechanism-based models, in which many 
parameters are often not identifiable based on data typically 
collected in clinical trials. The Bayesian approach is more 
commonly applied in toxicology research, in which the algo- 
rithms and software available to implement such approaches 
have been widely adopted. 1213 However, examples of more 
clinically oriented PK/PD models developed using Bayesian 
approaches have also been reported. 14-16 MCMC Bayesian 
algorithms offer distinct advantages over more commonly 
used deterministic estimation algorithms. It does not suffer 
from the notorious numerical instability that often occurs with 
linear approximation algorithms (e.g., FO and FOCE) and 
complex highly dimensional PK/PD models. 17 In addition, 
formal utilization of prior knowledge of PK/PD relationships 
using MCMC Bayesian algorithms allows leveraging of the 
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Figure 4 Clinical trial simulation of 4|3HC following 14 days of once-daily treatment with potential inducers of liver CYP3A4. (a) Mean levels 
of 4pHC vs. time, (b) Probability of detecting a difference relative to placebo vs. sample size on day 14. (c) Fold induction in CYP3A4 vs. the 
sample size to achieve 80% probability. 4pHC, 4p-hydroxycholesterol; CYP3A4, cytochrome P450 3A4. 



wealth of prior knowledge for model development. Thus, a 
model developed using this approach represents not only 
knowledge gained from the current analysis dataset but also 
potentially all prior knowledge. 

To be able to accurately estimate oral and i.v. MDZ expo- 
sures, separate compartments were needed to represent 
CYP3A4 interaction in the liver vs. the gut. Interestingly, the 
fold induction of CYP3A4 estimated for rifampin in the gut 
was much greater than that in the liver. This is consistent 
with reported effect of rifampin on verapamil, where oral 
clearance is elevated 32-fold and systemic clearance only 
increases 1.3-fold. 18 To accurately estimate the effect of 
ketoconazole on both MDZ and 4pHC, separate parameters 
for CYP3A4 K were necessary, with the K for 4pHC esti- 
mated to be 400-fold higher than that of MDZ. This finding 
could be explained by the potential of 4pHC to be generated 
through metabolism by other metabolic enzymes besides 
CYP3A4. This is consistent with the observation in a clini- 
cal study that propiverine treatment produces an increase in 
MDZ exposure, suggesting CYP3A4 inhibition, and a signifi- 
cant increase in 4pHC, suggesting induction of an alterna- 
tive pathway (non-CYP3A4/5) for formation of 4pHC. 19 The 
half-life of 4pHC was estimated to be 10 days, only slightly 
less than the half-life of 17 days commonly described in the 
literature. 310 However, the model-based estimation of the 
4pHC half-life is likely to be reliable for healthy subjects given 
the robustness of model-based approaches for estimating 



physiological rate constants compared with the noncom- 
partmental approaches typically used in analysis of clinical 
study data. There were differences in the MDZ PK param- 
eters between the posterior mean estimates and the priors, 
based on the meta-analyses by Yu etal. 20 and Li etal. 2 \ par- 
ticularly for V max and bioavailability. The meta-analyses by Yu 
et al. 20 and Li et a\ 2A were conducted using aggregate level 
(i.e., summary level) data, whereas the current analysis is 
conducted using individual-level data, and this may explain 
some of this difference. In addition, differences between the 
population sample used in the current analysis compared 
with the meta-analyses by Yu et al 20 and Li et a\ 2A may also 
be partly blamed. 

4pHC has been proposed as a potentially useful bio- 
marker for early detection of CYP3A4 inhibition and induc- 
tion in clinical development. The PK/PD model we have 
developed based on the clinical study in healthy subjects 
allows comparison of 4pHC vs. MDZ, the "gold standard" of 
CYP3A4 probe substrates, in realistic and clinically relevant 
scenarios. The simulations performed with the model sug- 
gest that dynamic range of 4pHC for detection of CYP3A4 
induction or inhibition is significantly narrower compared 
with that of MDZ. In particular, 4pHC is rather insensitive for 
the detection of CYP3A4 inhibitors; for example, an inhibi- 
tor that increases MDZ AUC by 20-fold would result in only 
a 20% decrease in 4pHC after 14 days of dosing. There- 
fore, 4pHC is unlikely to detect CYP3A4 inhibition for the 
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Figure 5 Clinical trial simulation of 4pHC following 14 days of once-daily treatment with potential inhibitors of CYP3A4. (a) Mean levels of 
4(3HC vs. time, (b) Probability of detecting a difference relative to placebo vs. sample size on day 14. (c) l/K. at C min vs. sample size for 80% 
probability, (d) Summary table of relevant parameters for CYP3A4 inhibitors. Drugs B and C have PK properties and CYP3A4 K. similar to 
ketoconazole. 4pHC, 4p-hydroxycholesterol; CYP3A4, cytochrome P450 3A4. 



majority of new chemical entities that are routinely devel- 
oped, unless the sample size is sufficiently large to detect 
small changes in 4pHC. 

In deciding to incorporate 4pHC as a biomarker for either 
inducers or inhibitors in clinical development, the sample 
size of the study must be taken into account. This is particu- 
larly true for inhibitors, due to the lower sensitivity of 4pHC 
for detection of CYP3A4 inhibition. We have used the PK/ 
PD model to simulate 4pHC in trials with different sample 
sizes, with both inducers and inhibitors of CYP3A4. For very 
strong inducers such as rifampin, there would be a high 
probability of detecting 4pHC elevation, even by comparing a 
single arm (N ~ 6) of a typical multiple ascending dose study. 
Detection of the signal could be enhanced by combining 
data from several panels of a multiple ascending dose study. 
For detecting the weak inducers (e.g., 1.5-fold induction), a 
much larger sample size (n > 60) would be needed to detect 
the 4pHC elevation with high probability. These simulations 
suggest that 4pHC could be a useful biomarker for detection 
of CYP3A4 induction for a new chemical entity early in the 
clinical development program. However, the utility of 4pHC 
to detect CYP3A4 inhibition is much more limited. For an 
inhibitor similar to ketoconazole, with an \I\IK. at C of 0.1 1 

' L J i mm 

(Figure 5), administered with a dosing regimen of 200 mg 
b.i.d., it would be feasible to expect detection of 4pHC reduc- 
tion in a typical multiple ascending dose study. Not sur- 
prisingly, the same inhibitor administered as a 400 mg q.d. 



regimen (Figure 5) would require a larger sample size. Most 
CYP3A4 inhibitors that are now in clinical development are 
likely to be weaker than ketoconazole. In addition, it is likely 
that the accuracy and precision of the 4pHC assay would 
make the identification of these relatively small changes dif- 
ficult to determine with any degree of confidence, making 
the potential for detecting false-negatives or false-positives 
unacceptably high. Therefore, the utility of 4pHC to be used 
for detection of CYP3A4 inhibition early in clinical develop- 
ment is low. 

There have been several clinical trials conducted to evalu- 
ate the clinical utility of 4pHC as a biomarker for inducers 
or inhibitors of CYP3A4. These trials have utilized commonly 
known perpetrators of CYP3A4 interactions; inducers such as 
rifampin, efavirenz, and carbamazepine; 6 ' 7 ' 9 or inhibitors such 
as ketoconazole, itraconazole, and ritonavir. 5 ' 922 Reports of 
the utility of 4pHC in the early detection of CYP3A4 induction 
or inhibition in clinical development have yet to be published. 
The simulations performed in this study demonstrate the lim- 
ited utility of 4pHC for detection of CYP3A4 inhibition in early 
clinical development. However, detection of CYP3A4 induc- 
tion appears to be possible with the sample sizes typically 
used in phase I clinical trials. In either case, incorporation of 
4pHC in phase I studies should be relatively straightforward, 
and significant changes detected in these studies could be 
used to trigger a more definitive MDZ interaction study later 
in clinical development. 
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METHODS 

Development of the PK/PD model for 4pHC 

PK and PD data from a clinical study of the effects of rifampin 
and ketoconazole on the metabolism of endogenous markers 
of CYP3A4 activity, and on MDZ PK, in 34 healthy subjects 
aged between 18 and 45 years, were used in the analysis. 23 
All subjects were randomized to receive oral (2mg) and i.v. 
(0.4 mg) MDZ on day 1 or 2, and on day 1 5 or 1 6 after receiv- 
ing q.d. rifampin 600 mg, ketoconazole 400 mg, or placebo for 
2 weeks. Subjects were followed up until day 30. Serial blood 
samples for MDZ PK were collected on days 1,2, 15, and 1 6. 
Predose blood samples were collected from day 1 to 30 for 
assessment of 4pHC, 4ocHC, and total cholesterol. The data 
used in the PK/PD analysis included 1 ,600 MDZ plasma con- 
centrations, 606 plasma total cholesterol concentrations, 606 
plasma 4pHC concentrations, and 605 plasma 4ocHC con- 
centrations. The lower limit of quantification was 0.3 nmol/l 
(0.1 ng/ml) for MDZ plasma concentration and 2ng/ml for 
total cholesterol (52 umol/l), 4pHC (5 nmol/l), and 4ocHC (5 
nmol/l). Observations below the lower limit of quantification 
were ignored in the analysis. 

Development of the nonlinear mixed effects model was con- 
ducted through a MCMC Bayesian estimation approach with 
both informative priors and noninformative priors (Table 1) 
using the NONMEM software system (version VII, level 2.0, 
ICON Development Solutions, Ellicott City, MD). Intersubject 
variability was assumed to be log normally distributed for all 
fixed effect parameters except for the bioavailability param- 
eters, which were assumed to follow a logit distribution con- 
strained between 0 and 1 . A normal distribution was assumed 
for the prior fixed effects and an Inverse-Wishart distribution 
for the prior intersubject variance. No prior distributions were 
used for the residual variances. The PK models and prior dis- 
tributions for fixed and random effects for ketoconazole and 
MDZ are based on previously published meta-analyses. 2021 
Key model equations have been included in Figure 1, and 
descriptions of model parameters can be found in Table 1 . All 
model equations and parameters can be found in the NON- 
MEM control file included in the Supplementary Data. 

Thirty thousand iterations of each MCMC chain were used 
during the burn-in phase and 1,000 iterations for the sta- 
tionary phase. Plasma ketoconazole concentration was not 
measured in the clinical study; therefore, highly informative 
priors for ketoconazole PK parameters were used (Table 1). 
The effect of CYP3A44 induction was incorporated as previ- 
ously described: as an exponential function that increases 
with time and would be expected to reach steady state by 
-20 days of treatment, based on half-life of CYP3A4 (70h). 11 
Therefore, the dose response of inducers is not incorporated 
in the model, but rather the fold induction can be used as 
input to predict the change in 4pHC as a function of time. 
Convergence of the model parameters and objective func- 
tion values was evaluated by visual inspection of the "cat- 
erpillar" plots for three independent MCMC chains and by 
examination of the Gelman and Rubin statistics. 24 Gelman 
and Rubin convergence statistics for the three MCMC chains 
were generated using the Bayesian Output Analysis Program 
(BOA) version 1.12 (ref. 25) in S-Plus (TIBCO Spotfire S+ 
8.2;TIBCO Software, Princeton, NJ). 



Evaluation of the PK/PD model for 4pHC 

The posterior predictive performance of the PK/PD model 
was evaluated by generating predictions for the 34 subjects 
included in the analysis using 500 of the posterior sets of 
population mean parameter estimates. The 500 sets of 
parameter estimates were randomly sampled from the 1 ,000 
sets generated during the stationary phase of the MCMC 
distribution. Predictions were generated using NONMEM. 
AUCs were calculated for the observed and predicted con- 
centrations using the trapezoidal method. Summary statis- 
tics and plots were produced using S-Plus (TIBCO Spotfire 
S+ 8.2; TIBCO Software). 

Simulation of dynamic range of 4pHC vs. MDZ 

Simulations to compare the dynamic range for 4pHC vs. 
MDZ were performed by assuming 14 days of treatment 
with either inhibitors or inducers of CYP3A4. Two hundred 
and fifty inhibitors were simulated with CYP3A4 K. values 
ranging from 1 x 10 -7 nmol/l to 10 umol/l. Inhibitors were 
administered as 400 mg q.d. regimens with population 
mean PK parameters identical to those for ketoconazole. 
Two hundred and fifty inducers were simulated that ranged 
in induction potential for CYP3A4 from 0- to 100-fold and 
were assumed to achieve maximal induction with the same 
rate constant as rifampin. It should be noted that a one- 
fold induction is equivalent to a 100% increase in CYP3A4 
enzyme activity. All simulations were conducted using R 
version 3.0 (http://www.r-project.org/). 

Clinical trial simulations of CYP3A4 induction and inhibition 

Clinical trials were simulated to evaluate the sample sizes 
needed to detect changes in 4pHC after 14 days of q.d. dos- 
ing of either CYP3A4 inhibitors or CYP3A4 inducers. The 
change in 4pHC was evaluated relative to a placebo group of 
equivalent sample size. Five different CYP3A4 inhibitors with 
different [l\/K. at C min were simulated with sample sizes from 
6 to 100. To evaluate the effect of [f\/K. at C min on the sample 
size for 80% power, different [I\IK. properties were created by 
varying the dose and regimen of the five different CYP3A4 
inhibitors. Trial simulations for inducers were conducted in a 
similar manner, with inducers that increase CYP3A4 activity 
by 1 .5- to 4.6-fold. The probability of successfully detecting a 
significant change in 4pHC was defined as the percentage of 
500 simulated trials in which the 95% confidence interval for 
the treatment group does not overlap with 95% confidence 
interval in the placebo group. Type 1 error was estimated by 
evaluating the probability that the 95% confidence intervals 
between two groups treated with placebo do not overlap with 
sample sizes of eight subjects. All simulations, summary sta- 
tistics, and plots were generated using R version 3.0 (http:// 
www.r-project.org/). 
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Study Highlights 



WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

y The results of clinical studies of 4(3HC as a 
biomarker of CYP3A4 drug interactions are 
consistent with simulations presented in this 
article. Few publications have demonstrated 
the development or application of Bayesian PK/ 
PD modeling. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

y (i) How to use a Bayesian nonlinear mixed ef- 
fects approach to supplement gaps in PK and/ 
or PD data when developing a mechanism- 
based model? (ii) What is the potential utility of 
4(3HC as a biomarker for detection of CYP3A4 
drug interactions in clinical studies? 

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

1/ (i) A Bayesian approach to PK/PD modeling 
can effectively leverage prior knowledge to 
supplement PK/PD gaps in the analysis data- 
set, (ii) 4(3HC may be a useful biomarker for 
moderate to strong inducers, and strong inhibi- 
tors, of CYP3A4 in typical clinical studies. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

y Bayesian approaches will improve the man- 
ner in which prior information is used in PK/ 
PD model development. This article provides 
guidance for the incorporation of 4(3HC as a 
biomarker of CYP3A4 interactions in clinical 
development. 
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